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Abstract 

Drying a droplet of colloidal dispersion can result in complex pattern 
formation due to both development and deformation of a skin at the drop 
surface. The present study focus on the drying process of droplets of 
colloidal dispersions in a confined geometry where direct observations of 
the skin thickness are allowed. During the drying, a buckling process is 
followed by a single depression growth inside the drop. The deformation 
of the droplet is found to be generic and is studied for various colloidal 
dispersions. The final shape can be partly explained by simple energy 
analysis based on the competition between bending and stretching defor- 
mations. Particularly, the final shape enables us to determine precisely 
the critical thickness of the shell for buckling. This study allow us to val- 
idate theory in 2D droplets and apply it to the case of 3D droplets where 
the thickness is not accessible by direct observation. 

1 Introduction 

Drying droplets of complex fluids as polymer solutions or colloidal dispersions 
involves a large number of microscopic phenomena: solvent diffusion, transfers 
at the vapor/medium interface, skin formation then skin deformation. Over the 
last decade, there has been a great deal of scientific and technological interests 
in studying the flow and deposition of materials in drying droplets m. As 
an example, the fabrication of dried milk requires the transformation of liq- 
uid droplets into a powder form using nozzles; also liquid food concentrate is 
atomized into droplets using spray-drying processes [24l |l^ . Particularly 
functional properties have to be controlled such as dispersibility, and solubility 
that are essential for the capability of the powder to be re-hydrated. There- 
fore, the droplet morphology, or the droplet size changes have to be controlled 
and the effects of various initial or environmental conditions on a single droplet 
development have to be understood. 

The concentration at the free surface induced by desiccation can cause skin 
formation of the film [TZ and is often responsible for strong distortions of the 
droplet [m [21] . Deformation of nearly spheroidal droplets were considered in 
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the case of polymer solutions [TF or colloidal dispersions [191 [29]: the various 
shapes have been related to buckling process of a skin formed at the drop sur- 
face. Such deformations are not specifically related to drying droplets and can 
be observed at different scales: capsules formation |26] or lock particles for- 
mation when controlled shell buckling that forms spherical cavities [2 7 . However 
in spheroidal structure the thickness of the skin is not directly accessible gen- 
erally. Thus, we focus here on a single drop of a colloidal dispersion confined 
between two circular glass plates and left to evaporate (see figure [T]). In this con- 
figuration, diffusion processes and drying kinetics have been carefully studied 
[3 [161 [m M [3 [10] • Moreover this geometry allows us a direct observation and 
quantification of the skin growth near the drop/ air interface by imaging analysis. 
When its structural strength is sufficient, the skin can withstand internal stress 
and can be regarded as an elastic shell. By appropriate selection of physical 
quantities, constrained shrinkage can result in buckling instability, causing the 
circular shell to inversion of curvature as shown in sequence in figure [2^: the 
resulting single depression grows and is continued by an invagination tube pen- 
etrating inside the drop. The complex pattern formation is studied for various 
colloidal dispersions and drop sizes. The generic shape is explained by simple 
energy analysis. Particularly, the critical thickness of the shell for buckling is 
precisely determined by the final shape of the drop. In addition, this study al- 
low us to validate theory in 2D droplets and apply it to the case of 3D droplets 
where the thickness is not directly accessible. The case of three-dimensional 
droplets and large Peclet numbers corresponds to common industrial processes 
(spray-drying consisting in rapidly drying aerosols to manufacture dry powders). 

2 Experimental 

Our experiments were performed with various stable aqueous dispersions (see 
Table [T]). (i) silica particles Ludox HS-40 (solid volume fraction <p = 0.22, parti- 
cle diameter 2a = 15 ±2 nm, polydispersity index ~ 0.16) and diluted dispersion 
of Ludox TM-50 (solid volume fraction cf) = 0.12, particle diameter 2a = 22 ± 2 
nm, polydispersity index ~ 0.2 [9 ) from Sigma- Aldrich; (ii) Nanolatex particles, 
polystyrene (solid volume fraction (/) = 0.25, particle diameter 2a = 25 nm, glass 
transition temperature = 100°C) provided by Rhodia Recherche (Aubervilliers, 
France); (iii) Green fiuorescent silica particles sicastar-greenF (solid volume frac- 
tion (j) = 0.10, particle diameter 2a = 50 =b 3 nm, polydispersity index ~ 0.2, 
absorbance: Xexcitation = 475 nm and emission: )^emission — 510 nm) commer- 
cially available from Micromod. Most of the experiments were performed with 
the last system. The colloidal particles are adequately polydispersed to avoid 
cristallization. In addition during the drying process, the particles density in- 
creases resulting in the formation of a gel phase, defined as a porous matrix sat- 
urated with water. Dilution of dispersion was possibly done by adding deionized 
water (quality Milli-p). The confined geometry consists of a thin cell made of 
two circular, parallel and horizontal glass slides of radius Rs (figure [3]). It has 
been shown that the drying kinetics can be controlled by the slide size and is 
not strongly affected by the relative humidity of the surrounding [3; experiments 
were conducted at room temperature and relative humidity (~ 40%). The slides 
are lubricated with a thin film of silicon oil VIOOO (viscosity 1 Pa.s). This treat- 
ment prevents from the pinning process of the contact line [21^. A drop of the 
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Figure 1: (a) Set-up: a droplet of solution is sandwiched between two circular 
glass slides in side view (left); the cell is illuminated by transmitted light, (b) 
Sketch in top view of a quarter of cell showing the formation of a densely-packed 
particles at the liquid- vapour interface during solvent removal. 

dispersion is placed on the bottom substrate with a micropipet. Then the upper 
glass slide is carefully placed to squeeze the drop (see figured^). The gap (5, 
between the sHdes is controlled using three thin spacers of controlled thickness. 
In all experiments the radius of the glass slides is kept to = 8 mm and the 
gap is constant, 6 = 100 ± 5/im. During the drying process, images of the drop 
evolution are recorded at different times and show a dark ring corresponding to 
the meniscus at its periphery as shown in figure [2^. The high contrast of the 
images allows an accurate detection of the drop edge and provides variations 
with time of both the periphery length and surface area of the drop. Fluorescent 
microscopy was performed using a DM2500 Leica microscope with objective 5x 
and 50 X magnification for global view (drops deformation) and close-up view 
(envelope thickness measurements) respectively. The measurement of the size of 
the densely-packed region at the drop periphery, envelope thickness, have been 
investigated by fiuorescence microscopy (see figured]). Focusing on the air/drop 
interface, in the middle of the cell (z = | in figure |4^), the size of the envelope 
thickness, /i, can be estimated from the fiuorescence intensity profile with an 
accuracy of 5/im (figure IDd). 

3 Results 

The deformation of the drop during the drying process exhibits different con- 
figurations shown in the sequence in figure [2^ and in the superposition of the 
drop profiles in figure [2]3. The successive configurations adopted by the drop are 
driven by the volume decrease that is induced by solvent removal. In a first stage 
the drop progressively shrinks with solvent evaporation (figure [2^(1), (2)). Then, 
contrasting with the case of a pure solvent, the droplet stops shrinking isotrop- 
ically: a sudden inversion of curvature occurs (figure [2^(3)). Consequently, the 
depression is continued by an invagination tip that deepens with time (figure 
[2^(4), (5)) until the formation of a bean-shaped material that cracks in figure 
[2^(6)). These configurations are detailed in the following. Most of quantities 
used are shown in figure [3l 

3.1 Isotropic shrinkage and shell formation 

During the drying process, the gas that escapes from the liquid/air interface is 
driven outward from the cell by diffusion [5]. Direct observation of the drop 
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evolution with time shows a progressive and isotropic shrinkage in a first stage 
(figure [2^ (from (1) to (2)), figure [^d): this first stage is named as configuration 
1. As a result both periphery length and surface area decrease steadily with 
time [21] (figure [2]i). In the following we suppose that the radius decreases 
Hnearly with time as: R{t) = Rq — Jt where Rq is the initial drop radius and J 
is the evaporation rate. In addition time variations of dimensionless periphery 
length are shown in figure [2^ for different drops. In particular, the duration of 
the shrinkage process increases when the initial volume fraction of dispersions 
dereases (this is observed in the case of pure and diluted Ludox HS-40). 

Starting from suspended compounds uniformly distributed within the droplet, 
particles are advected to the drop periphery due to solvent removal (figure [TJd). 
The radial flow, induced by the evaporation rate, gives the typical velocity scale 
characterizing the transport of particles. A Peclet number can be defined as the 
diffusion time to = ^ (characteristic time for the diffusion of a particle along 

a distance Rq) divided by the evaporation time t^; = as follows: Pe = ^^^^ 
with the diffusion coefficient Dn = = 8.5 x 10~^^m^.5~^ using the Stoke 
Einstein relation {rj is the viscosity of the suspending fluid) and the evaporation 
deduced from the time variation of the periphery length, J ~ 10~^m.s~^ . In 
our experimental conditions: 10 < Pe < 10^, diffusion is slow compared with 
the rate of evaporation. As a result, the transport of particles to the drop- air 
interface strongly suggests the formation of a porous envelope during the evapo- 
ration process [13^ i28^ . This envelope thickens with time as shown by the direct 
measurements in figure |4]i. 

Due to the drop shrinkage, the structure of the envelope is possibly modified 
by a redistribution of particles. In this hypothesis, the particle volume fraction 
in the envelope increases until a maximum value close to the random close 
packing (j)c- Applying mass conservation of particles the relative thickness /i, of 
the shell of radius R can be determined during isotropic shrinkage [29^. 

Figure [1] allows us to plot the time variation of the envelope thickness from 
the time variation of the drop radius. Measurements of the envelope thickness 
are well fitted by this expression in the isotropic regime (figure |3]i). Due to 
the meniscus the thickness is possibly underestimated which could explain the 
discrepancy with expresion[Tl The thickening of the envelope does not impede 
solvent loss since the transport takes place through the envelope in accordance 
with the Darcy law: 

T] h 

where AP is the pressure drop across the shell thickness /i, required to 
produce a given fluid flow and k is the permeability of the porous shell estimated 
using the Carman-Kozeny relation [4 . 

When the structural strength of the envelope is sufficient, the envelope can 
withstand internal stress. At this time, particles do not redistribute in the 
envelope that will be considered as an elastic shell. The material is assumed to 
be homogeneous and characterized by an elastic modulus E] E is supposed to 
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Figure 2: (a) Digitized images taken during a typical deformation of a drop of 
colloidal dispersion (Ludox HS-40 (diluted), Rq = 2.9mm, Rs = 8mm). (b) 
Time variation of the relative volume of the drop (measurements); at the buck- 
ling onset a indentation depth, can be defined and characterize the growth of a 
depression, (c) Superposition of drop profiles during the drying process: dura- 
tion elapsed between two images is 1 minute, (d) Time variations of periphery 
length (crosses) and surface area (circles) of the drop; the inversion of curvature 
takes place at time t^, says the buckling time. The dashed line corresponds to 
the linear decrease of the drop periphery length with time, (e) Time variations 
of the dimensionless periphery length for various drops. 
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Figure 3: Mechanical stress field in cylindrical coordinates (left part of the 
sketch) and definition of physical quantities in the drying droplet (right part of 
the sketch). 



be constant during the deformation since the gel is continuously drained by the 
solvent. 

In addition the mechanical stress is assumed to be homogeneously distributed 
in the elastic shell. Observing that the shell thickness is much smaller than the 
other two spatial lengths (drop radius R and cell gap S), the following assump- 
tions can be made. Let us consider the three components of the mechanical 
stress field in the elastic shell, in cylindrical coordinates: the ortho-radial, radial 
and normal components of the stress field are denoted by aoo^ cFrr and cFzz re- 
spectively (see figure [3]). Equilibrium of the corresponding internal forces in the 
elastic shell results in the following relations: ^ >> 1 and § >> 1 

since ^ << 1 during the elastic deformation of the envelope. Therefore the 
main component of the stress in the shell is the ortho-radial one. This quan- 
tity is responsible for the buckling process of the shell. Indeed, above a critical 
thickness, deformations occur mainly by bending, which is much less energetic 
than stretching. Therefore a buckling process occurs leading to an inversion 
of curvature of the shell (figure [^i and figure [Sb). In the following suffix B is 
related to quantities taken at time ts- Expression [2] gives the relation between 
the shell thickness hs^ and the pressure drop APb across /ib at the onset of 
buckling: 

hB = -^- (3) 

Firstly expression [T] at t = ts gives as a function of Rq and hs- Replacing 
hshy the expression[3l ^ can be expressed as a function of Rq and the pressure 
drop across the shell at the onset of buckling, APg. In this way APb is a 
parameter and can be estimated for each system studied. ^ slightly increases 
with the initial drop radius as shown by the measurements for different colloidal 
systems (figure [5ti). These results are close to the linear increase of Rb with Rq 
obtained in spherical drops by Tsapis et al [29] . Secondly, combining expressions 
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(a) top view (b) 
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Figure 4: (a) Fluorescence images in top view (dispersion of Green fluorescent 
silica particles), (b) Corresponding intensity profiles across the air /drop in- 
terface, along the dashed line in (a); in the referential bound to the drop/air 
interface, the increase of the envelope thickness h{t), can be measured with time, 
(c) Menisci at the drop edge in side view z = and z = S denote positions with 
the substrates {5 = 100 ± 5/im); the contact angle of the aqueous dispersion on 
the surface is close to 60^. (d) Ratio of the envelope thickness /i, to the drop 
radius as a function of time: crosses are measurements deduced from the 
intensity profiles in (b) {Rq = 2.9mm, Rs = 8mm), dashed line is theory in 
accordance with expression [1] in the isotropic regime. 

[T] and [3] the critical thickness at the point of buckling is plotted in figure [Sb; 
measurements of the thickness using the method described in figured! The order 
of magnitude of APb is found to be in agreement with the values obtained by 
Tsapis et al.[29\. Values are listed in Tabled] for different systems. 

3.2 Depression growth 

As shown previously, the shell stops shrinking {R = Rb) and a buckling process 
results in the formation of a depression in a region of the shell (figure and 
figure [6^). In the following, the existence of the depression in the shell is referred 
as configuration II. A part of the elastic energy is released by the formation 



Table 1: Main characteristics for the samples considered in the experiments: 
particle diameter 2a (and polydispersity is known), particle volume fraction (/), 
permeability of the gel /c, and critical pressure drop across the shell at the onset 
of buckling APe, deduced from the measurements. 



2a{nm) (j) fe(xlO~^^m^) APBjkPa) 



Ludox HS-40 


15 ±2 


0.22 


1.2 


9 


Ludox HS-40 (diluted) 


15 ±2 


0.10 


1.2 


9 


Ludox TM-50 (diluted) 


22 ±2 


0.12 


2.4 


5 


nanolatex 


25 


0.25 


4.4 


3 


Green fluo silica 


50 ±3 


0.10 


12 


0.8 
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Figure 5: (a) Ratio of the radius of the sheh at the onset of buckhng to the 
initial drop radius as a function of Rq. (b) Ratio of the critical envelope 

thickness to the drop radius at the point of buckling as a function of Rq. 
Dots are measurements and lines are fits to the theoretical expressions [T] and [3l 

of this depression (figure [6^1). The depression growth is characterized by the 
increase of the indentation length e, and happens in a region of linear size £. 
Measurements of ^ as a function of e are shown in figure [GJd. These results 
are well fitted by the geometric relation: £ = 2y/2RBe — approximated by 
2yj2RBe to first order in e, as a consequence of the inversion of curvature. 
The elastic energy of configuration II is the sum of three contributions: 

Ull = /7o + Ufold + Uinversion (4) 

1. The term Uq is the contribution of the elastic energy due to the part of 
the shell not inverted: it is considered as a reference energy. 

2. During the depression growth the elastic energy contains a main contri- 
bution coming from the creation of the folds A and B (see figure [6^): a 
fold is limited by asymptotes (asymptotes A and A' for A) and extends 
over a typical distance d. Assuming h/R « 1, which is actually realized 
during the drop deformation, the Foppl-von Karman theory for elastic 
shells [8j HSl H] allows us to obtain the energy due to the fold Ujoid- 

Minimizing the total elastic energy (bending energy and in-plane elastic 
energy) with respect to the distance d leads to d = {HRbY^'^ (this was 
described in references [151 HI] ) • Consequently the relevant curvature ra- 
dius is tana where a is the angle between the asymptotes limiting the 
fold (figure [6^). Therefore Ufoid ^ tana/2^, where s = 2d3 is the 
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Figure 6: (a) Top: superposition of digitized profiles showing the depression 
growth in a part of the shell; the whites dots locate the regions of highest 
curvature. Down: sketch of the depression formed by inversion of the circular 
cap: the indentation length e, happens in a region of size between A and B. 
The angle a between the asymptotes A and A' limit the fold; the quantity d is 
the lateral extension of the fold limiting the circular shell and the inverted cap. 
(b) Ratio of the lateral size of the depression, iexp^ to the radius of the drop 
at the onset of buckling Rb^ as a function of the indentation length e, to the 
critical thickness /i^, at the onset of buckling; £^ = 2^/2RBe — e^, £ = 2^2RBe. 
Hb — 12 ± 5/im. Above a critical indentation length, the quantity i becomes 
constant: £ = £c- (c) Sketch of the shell forming at the drop periphery: the 
geometrical quantities of the shell are defined taking into account the shell 
thickness h. 
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area of the fold[22j [53]. As the depression grows, the angle a increases, 
in accordance with tana/2^ ~ so does the energy Ufoid- Ufoid is 
due to both the increase of the indentation length and the increase of the 
shell thickness. Consequently, at each time denoted by t^, an indentation 
length e^, and a shell thickness /i^, can be defined. 

Ufoidihi, Ci) = 2co5Eh^/'^^j^ei (5) 

where cq is a parameter that only depends on the Poisson ratio of the 
shell. 

3. The term Uinversion corresponds to the change in elastic energy due the 
inversion of the circular section AB in accordance with figure [6t. In 
this way the arc of external radius Rb + hB^^ becomes compressed after 
inversion; also the arc of internal radius Rb — hB^^ becomes stretched 
after inversion. The length change makes this energy contribution linear 
in the deformation [23] . 

TT (h \ — T^^h^ -^/^ ^f{^ 

^ inversion yi^i-) ) — ^ Jz/Orij^ 3/2^* 

R^ 

where ci is a parameter that only depends on the Poisson ratio of the shell. 
Note that Uinversion{hi^ Ci) bccomcs negligible compared with Ufoid{hi, e^) 
for Ci » h. 

The energy (Uu — Uq) / {EK^S) is plotted as a function of the relative volume 



variation, 



Ay 



~ ^'^^^2^ — - — in figure [T] This illustrates the tendency of 
the elastic energy to increase during the shell deformation. 

Since the energy concentrated in the folds becomes enormous during the 
depression growth, another configuration, denoted in the following by config- 
uration III, is energetically more favourable. This is realized by assuming the 
angle a constant {a defined in figure [6^^). Hence, the lateral size of the depres- 
sion ^, stops increasing and becomes constant as shown by our measurements 
in figure [6)3 (region named III). Now, the elastic energy due to the fold only 
increases through the growth of the shell thickness. 

However, the decrease of the inner volume of the shell need another mode of 
deformation of the shell. Assuming a constant and keeping i constant, as shown 
by empirical observations, require an increase of the indentation length e (figure 
[6^1) and results in a length change of the inverted part of the shell between A 
and B. Let us consider the stretching energy due to the length change of the 
inverted part of the shell. Usually, when an elastic material is deformed the 
work is stored in the body as a strain energy. If the elastic deformation 
(strain), occurs in an elementary volume dV ~ dhi {£ being the length scale of 
the stretched region), the elastic energy expresses as: 

As a result, the energy of the inverted part due to the relative change in 
length u expresses as: 
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Figure 7: Dimensionless elastic energy as a function of the relative volume 
variation | \ (the lowest-energy state is plotted as a solid line and the higher- 
energy one is plotted as a dashed line); the plotrange starts at the buckhng 

process | |= At a threshold value, | |, configuration II becomes 
energetically more favourable than configuration III. The high stretching process 
occuring during the configuration III possibly causes weakening and breakage 
of the inverted part of the shell; in a realistic point of view the vertical dashed 
line is not reached during the deformation of the elastic shell. 



Ustretch(hi, Ci) = C2SEhiiu{ei)'^ (7) 

where C2 is a parameter that only depends on the Poisson ratio of the shell. 
The relative length change u{ei) differs from zero only for > Cc. Then, the 
length change is related to the indentation length b}0 u{ei) ~ ^ ^r^* 

The contribution U stretch governs the elastic energy of the configuration III 
as follows: 

Ulll = /7o + Ufold{hi, Cc) + Uinversion{hi, 6^) + Ustretch{hi, 6^) (8) 

Here, Ufoid{hi^ec) only varies by the thickness variation because i is con- 
stant. 

The energy {Uj 1 1 — Uq) / {Eh'^ S) is plotted as a function of the relative volume 
variation | ^ | in figure [71 for small deformations, Uu < Um and configu- 
ration II is energetically more favourable than configuration III. Note that this 
new configuration costs rapidly a high energy. 

Let us now predict the critical quantity e = ec, and consequently i = ic^ 
at which the transition between configuration II and III arises. Starting with a 
deformed shell characterized by an indentation depth e (configuration II in the 

■"^ Starting from a circular arc AB, the simplest perturbation corresponding to an elonga- 
tional shape is a part of an ellipse. We use the approximation by Ramanujan 3 for the 
circumference of an ellipse gives: 5) ~ 7r(3(a + 6) — ^y(3a^^^~S)(a+~3E)) where a and b are 
one-half of the ellipse's major and minor axes respectively. Taking a = Rb + and h = Rb-, 
the relative length change in length between such an ellipse and a circle of radius Rb expresses 

as: (^b+^z^^b) (Rb^^b) ^ to first order in e^. The relative length change related to 

the transition from configuration II to configuration III in [71 is then . 



11 



I 



I 



I 




f 



I 



Figure 8: Ideal scheme of the transition between different configurations of 
the elastic shell. Configuration II corresponds to the inversion of a part of 
the elastic shell characterized by an indentation length e and a lateral size i. 
Configuration II can be followed by two possible configurations: (i) configuration 
ir corresponds to the increase of both the lateral size and the indentation length; 
(ii) configuration III where the lateral size is the same than in configuration II 
and the indentation length is the same than in configuration IF; consequently 
a stretching process of the inverted part is needed. 

sketch figure [8|), an increase of the indentation length can result in two possible 
configurations: 

• both the indentation depth e and the lateral size i increase in accordance 
with configuration IF in figure [8|); 

• only the indentation length increases at fixed lateral size ^ = £c in ac- 
cordance with configuration III in [51 in this case a length change of the 
inverted region is needed. 

At the threshold, the energy of configuration IF has to equal the energy 
of configuration III. Assuming that the elastic energy due to the inversion, 
inversions IS dominated both by the energy due to the folds and the one due to 
the stretching process, we have: Ufoid{h,e-\-ec) = Ufoid{h,e)^Ustretch{h,e^ec) 
using the notations in figure [H Starting from configuration II with e « ec, 
it comes: Ufoid{h^ec) ~ Ustretch{h^ec) (at the threshold). Using relations [5] 
and [71 it comes directly that Cc is proportional to the thickness at the onset of 
buckling, hs, with a prefactor only depending on the Poisson ratio: Cc = S^hs- 
Taking into account numerical values of the parameters cq and C2 with a Poisson 
ratio equal 0.3 we obtain [2 : Cc ^ ^hs- In addition, since the transition to 
lower-energy results in a constant lateral size this length can be expressed 
as a function of the characteristics of the shell at the buckling onset as: ^ 



As a consequence the thickness of the shell at the onset of buckling hs^ can be 
precisely deduced from the measurement of the distorted drop at a macroscopic 
scale since: 



llB 




(9) 
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Figure 9: Measurements of the quantity ^ 7^ ^ range of initial drop 

radii i^o, and for different systems. 

Figure [9] shows the dimensionless quantity J^Rb ^ ^ ^ range of initial 
drop radii. Evaluation oilfiB using the relation [5] is related to the measurements 
of macroscopic quantities (/c, Rb) and is consequently more accurate than the 
direct measurement of hs- Typically, the error by direct measurement of hs 
is 20% while it is equal to 5% in the measurement of 32^- Note that the 
discrepancy between theory and experiment could be due to the systematic 
discrepancy between measurements of the enveloppe thickness and prediction 
as shown in figure 

The continuation of the deformation is shown in figure [TOk : the depression 
deepens and forms an invaginating tube (configuration IV in figure [TOb). As 
shown in figure [71 the stretching process of the shell is the main contribution to 
the elastic energy of configuration III. However the energy required to proceed 
with this deformation is rapidly enormous. Therefore a weakening or a breakage 
of the stretched part of the shell is suspected. 

Indeed, a relative length change larger than 20% as observed experimentally 
is often not feasible in such colloidal gels (figure [10]) . Let us compare the tensile 
stress in the stretched part of the shell with the cohesive stress in the shell. The 
critical stress dc, needed to separate two particles of radius a can be estimated 
from the van der Waals interaction in terms of equilibrium separation distance 
Zo, and Hamaker constant A, whose value depends on the surface chemistry of 
the particles[12 : (Jc ^ x a~^. 

This critical stress is easily balanced by the tensile stress in the inverted 
region of the shell during the relative length change u: Eu ^ cFc. Taking orders 
of magnitude for the Young modulus, E ~ lO^Pa, obtained by indentation 
testing, and A = 0.83 x 10"^°J, Zq ^ 0.2nm for Ludox HS-40, the relative 
length change needed to separate particles is found to be: = 0.18. This value 
corresponds to a relative volume variation of | AV/Vq |~ 0.52 (note that the 
same order of magnitude is obtained for the different colloidal systems studied). 
Therefore the strong stretching of the shell possibly causes rapidly decohesion 
of the gel as the indentation depth grows. Then, breakage or weakening of 
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Figure 10: The indentation depth e at fixed lateral size £ = £c turns into an 
invagination tip (superposition of profiles in (a)), (b) In configuration II, the 
length of the inverted cap, C, increases with the indentation depth e; the black 
line is the theoretical arc length of inverted cap of radius Rb indented by e (arc 
length = Rbcos~^{1 — Ae/RB))- In addition the lateral size £ increases with the 
indentation depth e. The vertical line delimits configuration II and configuration 
III. In configuration III, the lateral size keeps a constant value, £ = £c^ while 
the length of the inverted part, C, still increases. When the length change of 
the inverted part reaches a critical value, possible decohesion in the inverted 
part occurs leading to a transition to configuration IV. This last configuration 
is characterized by the progression of an invagination tube that deepens in the 
droplet. 



the inverted region is possibly followed by a healing process of the shell due to 
particles accumulation during drying. 

3.3 Case of spheroidal shells 

The case of three-dimensional droplets and large Peclet numbers corresponds to 
common industrial processes. As in the confined geometry (2D case), drying a 
droplet of a colloidal dispersion leads to strong distortions (figure [TTb). However 
in 3D case, measurement of the shell thickness is not possible by direct obser- 
vation. The previous study can easily be applied to three-dimensional droplets 
to obtain shell thickness at the buckling instability. 

The experiments were carried out with droplets of Nanolatex dispersion. The 
nearly spheroidal geometry is obtained by depositing droplets on a hydrophobic 
substrate; the surface is sprayed with Lycopodium spores which form a rough 
texture in which air remains trapped when a drop is deposited. In the first stage 
the droplet shrinks as described in the superposition of profiles in figure [TTb . 
Then an inverted region grows at the apex of the droplet (figure [TTb. b)[19j. As 
in the 2D case, the depression is continued by an invagination process. Similar 
arguments can be applied to the case of spheroidal drops as explained in refer- 
ences [22l [23] . In particular the energy due to the fold and the one due to the 
stretching to the shell express now as: 

U^Uh^.ei)=''^E^er (10) 

C/i',?,,,,(/.„e,) = #^^ef (11) 
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Figure 11: Droplet of a Nanolatex dispersion deposited on a super- hydrophobic 
surface, (a) Distorted droplet: digitized image showing a depression of lateral 
size Ic at the top of the droplet, (b) Superposition of dimensionless profiles 
(apex height vs. radius R) measured at different times by lateral imaging: 
the duration between two consecutive profiles is 300 s. 



Cq and C2 only depend on the Poisson ratio. From [TO] and [TT] it comes 
that Cc is proportional to the thickness at the onset of buckling Hb, such as: 

Since the geometric relation ic 2\flRB^c still holds, it comes: 

1 J 2/3 .2 1 72 

Consequently, using relation [121 nieasurements of the macroscopic quantities 
Ic and Rb allow us to evaluate with a good accuracy the critical thickness for 
buckling of spheroidal shells. In the case of a droplet of nanolatex dispersion 
(solid volume fraction 0.25, particle diameter 25 nm, glass transition tempera- 
ture = 100°C) dried at RU = 50% and T = 20°C, we find Hb = 25 ± 2fim. 



4 Conclusion 

Desiccation of drops of complex fluids can display large shape distortions related 
to the development of an elastic skin at the drop surface. We studied such 
mechanical instabilities in a confined geometry where direct measurements of the 
skin thickness are possible. The drop deformation appears to be generic and does 
not depend on the initial drop size: a depression due to an inversion of curvature 
in the elastic shell is formed and deepens in the drop. This deformation can 
be related to a mismatch between bending and stretching processes. From the 
comparison between the corresponding elastic energies, we obtain a relation 
between the critical thickness for buckling and the characteristic lengths of the 
final drop shape that are easily measurable. The validation of this method in the 
confined geometry (2D geometry) is then applied to spheroidal drops exhibiting 
identical distortions. 
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